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The experimental binding energies of single-particle and single-hole neutron states belonging to 
neutron shells that extend from N = 126 to 184 and 82 to 126 respectively, have been reproduced by 
solving the Schrodinger equation with a potential that has two components: the generalized Woods- 
Saxon (GWS) potential and the spin-orbit (SO) coupling term. The GWS potential contains the 
traditional WS potential plus a term (SU) whose intensity reaches a maximum in the nuclear surface. 
Our results indicate the existence of a explicit relationship between the strength of the SU potential 
and the orbital angular momentum quantum number i of the state. This dependence has been 
used to make reasonable predictions for the excitation energy centroids of states located inside and 
outside the neutron shells investigated. Comparisons are made with results reported in previous 
investigations. 

PACS numbers: 


1. INTRODUCTION 

Single-particle and single-hole neutron states have been previously investigated in the region around the doubly- 
magic 20S Pb nucleus [1-6]. The interaction of a neutron with the rest of the nucleus referred to as the core has been 
represented in these studies by a Hamiltonian containing a nuclear Woods-Saxon (WS) potential[7] and a spin-orbit 
(SO) coupling term. Although the same potential parameterization has been used in Refs. [1-4], a unique set of 
parameter values has not been found. In fact, none of the reported parameterizations, including the one published by 
Schwierz et al. [5] , has been able to produce an overall agreement between the predicted single-particle and single-hole 
neutron state energies and the corresponding experimental values. 

In this work, we have used the so called generalized Woods-Saxon (GWS) potential instead of the original WS 
potential]?], with the expectation of reproducing the experimental binding energies of single-particle and single-hole 
neutron orbitals that exist in the neutron shells N = 126 — 184 and 82 — 126, respectively. This potential contains 
the WS potential plus a term referred to as the surface (SU) potential that maximizes in the nuclear surface and is 
linearly proportional to the derivative of a WS function. It is a well known fact that the WS potential alone does 
not reproduce the energies of £ = 0 single-particle levels with enough accuracy when applied to a wide nuclidic region 
[8]. Using a formalism where the Schrodinger equation is expressed in terms of the Jacobi polynomials, Goniil et 
al. [8] have obtained the total potential influencing a single i = 0 neutron in a nucleus. The deduced potential (a 
GWS potential) contains a WS potential plus an additional term (a SU potential) that provides the flexibility to 
construct the surface structure of the related nucleus. The validity of the theoretical treatment presented in Ref. [8] 
is supported by the careful analysis of the analytical results obtained and the detailed discussion of the formalism 
prescriptions followed. 

The GWS potential has been the object of study in different works published in the literature. Hamzavi et al. [9] 
have obtained solutions of the Dirac equation for a nucleon experiencing the GWS potential under the relativistic 
spin symmetry limit. Bayrak et al. [10] have modelled the behavior of a neutral pion in a nucleus and obtained an 
analytical solution of the Klein-Gordon equation for a spin = 0 particle moving in the GWS potential field. Using 
the Dirac equation, Candemir et al. [11] have determined negative energy eigenvalues which satisfy the boundary 
condition for any k (spin-orbit operator eigenvalue) states of an anti-proton being subjected to the GWS potential 
in the pseudospin symmetry limit. Also, they have calculated the bound state energy eigenvalues of a proton moving 
under the GWS potential in the spin symmetry limit. These authors point out the usefulness of the SU potential to 
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examine the single-particle energy levels of a nucleon and an anti-nucleon since the interactions of the surface nucleon 
or anti-nucleon are very important in explaining the energy spectrum of the nuclei. 

The SU term of the GWS potential has been shown in some cases to generate an additional potential pocket in 
the nuclear surface region that is crucial to understand the elastic scattering of several nuclear reactions [12, 13]. In 
other cases, the SU potential does not induce this potential pocket as in the heavy ion elastic scattering analysis of I. 
Boztosun [14] where the inclusion, within the framework of the optical potential, of two small real SU potentials, was 
absolutely necessary to reproduce the elastic scattering data especially at large angles. In this case, although the SU 
potentials do not produce pockets, they generate interference effects that lead to an impressive reproduction of the 
experimental cross section oscillations [14]. The potential pocket mentioned here has also been obtained theoretically 
by Koura et al. [15] from the study of a refined WS potential that has a term that increases the freedom in the surface 
structure of the potential. 

C. Berkdemir et al. [16] used the Schrodinger equation with the GWS potential to generate the binding energies 
of several £ = 0 states of a hypothetical nucleus formed by a valence neutron plus an inert core of 56 nucleons. 
Unfortunately, the eigenvalue equation reported in this publication was later on shown to be physically incorrect [17]. 
However, the numerical calculations included in Ref. [17] show that the addition of the SU term to the WS potential 
allows the modification of the binding energies of £ = 0 states that otherwise would not be possible with the addition 
of a SO potential only. Recently, Bayrak et al. [18] examined the bound state solution of the radial Schrodinger 
equation with the GWS potential using the Pekeris approximation in terms of an appropriate boundary condition for 
arbitrary i states. From the closed form obtained for the energy eigenvalues, these authors calculated single-particle 
energies corresponding to a neutron orbiting around the 5e Fe nucleus that are comparable to the numerical results 
for 1 = 0. 

In this work, we use the GWS+SO potential mentioned above (GWS = WS + SU) to reproduce the experimental 
binding energies of zero and nonzero £ states in 209 Pb and 207 Pb. The WS and SO parameter values have been taken 
from the work published by Schwierz el al. [5]. For every state analyzed, we have obtained the strength of the SU 
potential that makes coincide the theoretical binding energy with the corresponding measured value. The experimental 
binding energies have been extracted from previously reported neutron transfer reaction experiments carried out on 
20S Pb [19-22] and the theoretical binding energies have been obtained by solving the Schrodinger equation with the 
computer program GAMOW [23]. The strength of the SU potential has been shown to be a function of the orbital 
angular momentum of the state. This dependence has been used to predict the excitation energies of some states that 
lie inside and outside the neutron shells investigated. 


2. THE TRADITIONAL WOODS-SAXON PLUS SPIN-ORBIT POTENTIAL 

2.1. Formalism 


As a starting point, we used a nuclear potential given by 


V(r) = V W s{r) + V S o(r), 

where the Vwsir ) and Vso( r ) potentials have the analytical expressions 

Vws 


Vws(r ) = - 


l + ex P (^ 


and 


t/ ( \ A \ 2 1 dV S o 
Vso(r) = --\ c 


j(j + l)-£{£+!)-- 


(1) 

( 2 ) 

(3) 


Here, A is a dimensionless parameter related to the depth of the SO potential and A c is the reduced Compton 
wavelength of the nucleon-core system given in fm units by 


Ac = — = 0.210019(l + -J^) , 

AV M core) 

(4) 

Vso(r) - Vs ° . , 

(5) 

1 +exp(-f-) 

R = r 0 A 1/3 , and R so = r S oA 1/ 3 . 

(6) 
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TABLE I: Experimental ( sp 'Bnij) and theoretical ( a ' p 'Bth,ntj) binding energies of single-particle neutron states in 209 Pb. The 
experimental binding energies were obtained from the expression s ' p ' B n ej = —3.9374 + E n ^(209), where the energy centroid of 
each orbital, E n ij( 209), is assumed to be equal to the excitation energy of only one excited state in 209 Pb since the spectroscopic 
factor ( C 2 S/(2j + 1)) measured by Kovar et al. [20] for each 2(19 Pb state, is close to 1. 


orbital 

. 1 

E n ij (209) 

C 2 S/(2j + 1) 

3p B nej 

sp B thMj 

(n<j) 


(MeV) 


(MeV) 

(MeV) 

2*79/2 

4 

0.000 

0.83 

-3.937 

-3.822 

lin/2 

6 

0.779 

0.86 

-3.158 

-2.699 

3^5/2 

2 

1.565 

0.98 

-2.372 

-1.915 

4Sl/2 

0 

2.033 

0.98 

-1.904 

-1.382 

2 <77/2 

4 

2.492 

1.05 

-1.445 

-1.117 

3^3/2 

2 

2.537 

1.07 

-1.400 

-0.945 


TABLE II: A comparison between the experimental ( s ' p 'B„ej ) and theoretical ( s ' p 'Bth,nej) binding energies of the lj’15/2 orbital 
in 209 Pb determined for different fragmentation schemes. The experimental energies were determined by evaluating Eq. (8). 
The theoretical energies were calculated with the computer program GAMOW [23] using the potential given in Eq. (1). The 
excitation energies of the 209 Pb fragments, E x (209), and their spectroscopic factors, C 2 S / (2j + 1), were extracted from 
an investigation of the 208 Pb(d,p) reaction at Ei a b(d) = 20 MeV [20]. For the scheme 1, the energy centroid E n tj(209) was 
assumed to be equal to the 200 Pb excitation energy E x ( 209)= 1.424 MeV. The schemes 2 and 3 include four 209 Pb fragments. 
The deuteron break up potential was used to generate the results corresponding to the scheme 3. 


Fragmentation E x ( 209) C 2 S/(2j + 1) E n tj( 209) 3 ' p B n ej a p 'B t h,nei 
Scheme (MeV) (MeV) (MeV) (MeV) 

1 

1.424 

1.00 

1.424 

-2.513 -2.369 

2 

1.424 

0.58 

1.782 

-2.155 -2.369 


3.052 

0.070 




3.556 

0.032 




3.716 

0.032 



3 

1.424 

0.77 

1.769 

-2.168 -2.369 


3.052 

0.09 




3.556 

0.04 




3.716 

0.04 




In these equations, /r represents the reduced mass of the neutron-core system, m n and M core are the neutron and 
core masses respectively, H is Planck’s constant divided by 2-7r, and c is the speed of light. The quantum numbers 
l and j correspond to the orbital and total angular momenta respectively of a single-particle or a single-hole state. 
Vws , R and a represent the depth, width and surface diffuseness respectively of the WS potential. The parameters 
corresponding to the SO potential are A Vso , Rso > and a so respectively. We selected the parameterization published 
by Schwierz el al. [5] 


V ws = Vo(l + j[-(N-Z + l) + 2]) , Vso = Vo , (7) 

and the parameter values: Vq = 52.06 MeV, k = 0.639, r 0 = 1.260 fm, rso = 1-16 fm, a = aso = 0.662 frn and 
A = 24.1. This WS + SO parameterization is applicable over the whole nuclear chart for nuclides between 16 O and the 
heaviest elements, and provides a relatively good description of the nuclear mean field leading to quality single-particle 
and single-hole spectra, nuclear radii, prediction of drip-lines, shell closures and other properties. 

To calculate the binding energy of single-particle neutron states, the 209 Pb nucleus (A = 209, Z = 82 and N = 127) 
is considered as a system formed by a valence neutron and a core of 208 nucleons, equivalent to the doubly-magic 
208pb nuc l eus . The reduced Compton wavelength (Eq. (4)) is calculated with m n = 1.0087 u and M core = 207.9767 u. 
For the determination of the single-hole neutron state energies, the 207 Pb nucleus is treated as a 20S Pb nucleus minus 
one neutron or, equivalently, a core of 207 nucleons ( M core = 206.9759 u) plus a neutron-hole. Since a single-hole state 
can be considered as an unoccupied single-particle state, the WS and SO widths and strengths, used to obtain the 
binding energies of hole (unoccupied) states in 207 Pb, are determined with the values A = 208, Z = 82 and N = 126. 
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TABLE III: Experimental ( 3 ' h 'B n ej) and theoretical ( s ' B t h, n ej) binding energies of single-hole neutron states in 201 Pb. The 
experimental binding energies were obtained from the expression s ' h 'B n ij = -7.3678 - E n ej( 207), where E n ej( 207) represents 
the excitation energy centroid of each orbital. The excitation energies of the 20 ' Pb fragments, E x ( 207), and their spectroscopic 
factors, C 2 S/(2j + 1), were extracted from Ref. [22] for the 3pi/ 2 , 2/ 5 / 2 and 3p3/ 2 states and from Ref. [21] for the rest of the 
states. 


orbital i E x ( 207) C 2 S/{2j + 1) 
(n£j) (MeV) 

E nlj { 207) 
(MeV) 

s.h. TD s.h. D 

-Esnij E>th,n£j 

(MeV) (MeV) 

3^1/2 

1 

0.000 

1.1 

0.000 

-7.368 -7.687 

2 / 5/2 

3 

0.570 

0.88 

0.570 

-7.938 -8.412 

3P3/2 

1 

0.900 

0.96 

0.900 

-8.268 -8.492 

1*13/2 

6 

1.629 

0.857 

1.869 

-9.237 -9.332 



5.990 

0.050 



2 / 7/2 

3 

2.334 

0.913 

2.913 

-10.281 -10.716 



4.572 

0.044 





4.765 

0.035 





6.370 

0.113 



I/ 19/2 

5 

3.415 

0.690 

4.014 

-11.382 -10.681 



3.660 

0.090 





5.410 

0.112 





5.620 

0.180 



17*11/2 

5 

7.010 

0.088 

8.077 

-15.445 -15.906 



7.590 

0.088 





7.960 

0.064 





8.220 

0.055 





8.540 

0.077 





9.220 

0.088 




The single-particle neutron orbitals studied in this work (2g 9 / 2: l*n/ 2 , lji5/2; 3d 5 / 2 , 45^, 2g 7 / 2 and 3 d 3 / 2 ) lie 
above the Fermi level, belong to a neutron shell that extends from N = 126 to 182, and all of them except the 
lji5/2 orbital have even £ values and positive parity. In contrast, the single-hole neutron orbitals (lh-n/ 2 ,l/ig/ 2 , 2 / 7 / 2 , 
l*i3/2> 3P3/2> 2/5 / 2 and 3 P1/2) are located below the Fermi level, all of them except the 1 hu / 2 level belong to the 
N = 82 — 126 shell, and all of them except the l*i 3 / 2 orbital have odd £ values and negative parity. We used the 
computer code GAMOW [23] with some minor modifications to calculate the binding energies of all these orbitals. 
The theoretical results obtained have been compared with the experimental single-particle and single-hole binding 
energies, s ' p B n y and s h B n £j, given respectively by 


s v B nlj = AM(209) - AM(208) - A M n + E nij ( 209) , (8) 

and 

s h B nij = AM(208) - AM(207) - A M n - E ntj ( 207) , (9) 

where AAf„, AM(209), AM(208) and AM(207) are the neutron, 209 Pb, 208 Pb and 207 Pb mass excesses equal to 8.0713, 
-17.6153, -21.7492 and -22.4527 MeV respectively. The energy E n (j( 209) (E n ^j (207)) represents the excitation energy 
centroid in 209 Pb ( 20( Pb) around which the strength of a single-particle (single-hole) neutron orbital (characterized 
by the quantum numbers n, £, and j) is fragmented. 


2.2. Results 

Table I depicts the theoretical binding energies, s p 'B t h,nij > determined with the GAMOW code [23] for the single¬ 
particle neutron orbitals 2 g 9 / 2 , 1 *n/ 2) 3d5/ 2 , 4si/ 2 , 2 g 7 / 2 and 2>d^/ 2 , using the potential given by Eq. (1). Although 
the lead isotopes have been extensively investigated since the beginning of nuclear physics, little experimental data 
on the fragmentation of these orbitals is available in the literature [19]. These orbitals have been predicted by the 
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shell model and observed in single-particle stripping reactions on 208 Pb [19]. Based on a meticulous investigation of 
the 208 Pb(d,p) 209 Pb reaction at Ei a b = 20 MeV carried out by Kovar et al. [20] and the results of different neutron 
stripping reaction studies compiled by Martin [19], the strength of each one of these single-particle states may be 
considered as mostly concentrated in one excited state of the 209 Pb nucleus. In order to determine each experimental 
binding energy s ' p B n cj shown in Table I, we assumed that each excitation centroid E n ^j( 209) used to evaluate Eq. 
(8), was equal to a particular 209 Pb excitation energy. 

Table II contains different values determined for the experimental binding energy of the l/ 15 / 2 orbital assuming 
three fragmentation schemes. From the study of the 208 Pb(a, 3 He) reaction at Ei a t,(a) = 183 MeV, Massolo et al. 
[21, 24] demonstrated the strong population of the lj’ 15/2 single-particle level at 1.424 MeV excitation energy in 209 Pb 
and the existence of fragments up to 5 MeV. In the first fragmentation scheme considered in Table II, the lj’ 15/2 
strength is assumed to be concentrated completely at 1.424 MeV excitation energy. The schemes 2 and 3 include 
three fragments, in addition to the 1.424 MeV state, that, according to Kovar et al. [20], exhaust most of the available 
I. 7 - 15/2 single-particle strength. Although these schemes include the same fragments, their corresponding spectroscopic 
factors, extracted from a DWBA analysis of the 208 Pb(d,p) reaction at Ei a b(a) = 20 MeV, are different because the 
deuteron break up potential was only included in the scheme 3. 

Table III displays experimental and theoretical binding energies of the single-hole neutron states in 207 Pb studied in 
the current work. Angular distributions of cross sections and analyzing powers corresponding to the ground (l/2 _ ), 
0.570 MeV (5/2“) and 0.900 MeV (3/2“) states in 207 Pb were measured by Matoba et al. [22] in a study of the 
208 Pb(p , d) 207 Pb reaction at Ei a b = 65 MeV. The spectroscopic factors obtained for these states indicate that the 
strengths of the 3p 1 / 2 , 2 f 5 / 2 and 3 p 3 / 2 orbitals are about 90 %. These spectroscopic factors are consistent with the pre¬ 
dictions of the particle-vibration coupling model published by Majumdar [25]. The experimental fragmentation data 
corresponding to the l*i 3 / 2 , 2 / 7 / 2 , I/ 19/2 and lhn/ 2 orbitals were extracted from a study of the 208 Pb( 3 He,a) 207 Pb 
reaction at Ei a b = 70 MeV carried out by Gales et al. [21]. 

The theoretical binding energies presented in Tables I and III for the 4 si/2, 3d 3 / 2 , 3d$/ 2 , 3pi/ 2 and 3^3/2 orbitals are 
consistent with those calculated by Wang et al. [6] using an elaborately constructed multi-step potential to approximate 
the WS + SO potential specified in Eq. (3). However, no agreement exists between any of the theoretical predictions 
and the corresponding experimental value, similarly to what has been reported in the literature [1—5]. As shown in 
Tables I and II, the theoretical binding energies of the single-particle neutron states are in general less negative than 
the corresponding experimental values. The opposite occurs for the single-hole neutron states displayed in Table III. 
This is a consequence of how Schwierz et al. [5] obtained the parameter values Vo = 52.06 MeV and n = 0.639, 
used in the current work to evaluate Eq. (7). These parameters were optimized by using a set of Vo and k values 
corresponding to single-particle and single-hole neutron and proton states present in the vicinity of the doubly-magic 
nuclei le O, 40 Ca, 48 Ca, 56 Ni , 100 Sn, 132 Sn and 208 Pb. 

3. THE GENERALIZED WOODS-SAXON PLUS SPIN-ORBIT POTENTIAL 

3.1. Formalism 

In order to correct the differences mentioned in Section 2.2, between experimental and theoretical energy values, we 
have replaced the traditional Woods-Saxon potential VVs(r) contained in Eq. (1) by the Generalized Woods-Saxon 
potential Vgws(t) [9-11, 16-18], given by 


Vgws{t) = V ws (r) + Vsu(r), (10) 

where Vsu( p ) is a surface potential, proportional to a derivative of a Woods-Saxon function, and given by 

( 11 ) 


Vsu(r) = -C I , ^^ J G, 


The parameter a was defined previously, R is given by Eq. (6) and the value of Ctj is calculated for each orbital when 
the theoretical binding energy coincides with the experimental value. 


3.2. Results 

Table IV contains the Cpj values obtained for all the single-particle and single-hole neutron states studied here. The 
positive (negative) signs of determined for single-particle (single-hole) states are, in part, a consequence of the 
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FIG. 1: A Ce = 1 / 2 ) - C u= e +1/2 ) vs. £{£ + 1). The solid line (Eq. (14)) was generated with the A Ce values of the s, d, 

g and i orbitals (circles). The A Ce values obtained for the p and f orbitals (solid triangles) lie close to a dashed line which is 
equal to the solid line minus 1.1142 MeV. The determination of the empty triangle was based on the excitation energy centroids 
measured by Gales [21] for the h orbitals. The square and the diamond are explained in the text. 


parameterization Eq. (7), and the associated parameter values [5]. We have also included in the table, the average 
Ce and the difference A Ce of every pair of Cpj values (j = t ± 1/2) corresponding to the same quantum number £. 
They are given by 


A Ce — C^j—e-i/ 2 ) — C^-e+1/2) , 


and 


Ce — [C(j—e~ 1 / 2 ) +C(j=e+ 172)]/2 


( 12 ) 


(13) 


Since the beginning of this investigation, we had the idea of finding a mathematical expression to connect the 
strength Cej of the SU potential with a physical quantity relevant to the problem. After a time consuming process of 
trial and error, we finally realized that when the differences A Ce (see Table IV) associated to £ = 0, 2, 4 and 6 (solid 
circles in Fig.l) are plotted as a function of the emblematic quantity £(£ + 1), they align along the straight solid line 
shown in Fig. 1. Surprisingly, the four solid circles displayed in Fig. 2 (corresponding to the Ce values for l = 0, 2, 
4 and 6 , plotted as a function of i(l + 1)) belong to the same parabola (solid curve shown in Fig. 2). These results 
together with the fact that Ce and A Ce are independent variables, support the existence of a formula that describes 
the strength Cej in terms of i{l + 1). The circles shown in Figs. 1 and 2, are nicely described by the solid curves 
given by 


A Ce = 0.1142£(£ + 1 ) - 0.2592 , (14) 

and 

Ce = 0.0058[£(£ + l )] 2 — 0.3851£(£ + 1 ) + 7.52 . (15) 

The squares, triangles and diamonds included in the figures will be explained in detail later. 

The impressive agreement observed between the circles and the solid curves shown in Figs. 1 and 2, allowed the 
determination of the Cej value corresponding to the lj’ 13/2 orbital. This state lies above the shell N = 126 - 184. 
Evaluating Eqs. (14) and (15) for l = 7, we obtained 7) =4 .080 MeV and /S.C^e— 7 ) = 6.127 MeV. In accordance 
with the notation used in Eqs. (12) and (13), C/j— 13 / 2 ) = Cu= 7 - 1 / 2 ) represents the Cp 3 value of the lj’ 13/2 orbital. In 
principle, C( ; ^i 5 / 2 ) may be any of the three Cej values displayed in Table IV for the IJ 15/2 orbital. It is important to 
realize that two C(j = i 3 / 2 ) values can be determined independently from Eqs. (12) and (13) when one of the possible 
C(j— 15 / 2 ) values is substituted into these equations. This independence allowed us to rule out the values Cq =15 / 2 ) 
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FIG. 2: C( = (Cj = i^\/2 + Cj=i- (-i/ 2 )/2 vs. £(l + 1). The solid line (Eq. (15)) was generated with the Ce values of the s, d, 
g and i orbitals (circles). The Ci values obtained for the p and f orbitals (solid triangles) lie close to a dashed curve which is 
equal to the solid line minus 8.9355 MeV. The determination of the empty triangle was based on the excitation energy centroids 
measured by Gales [21] for the h orbitals. The square is explained in the text. 


= -1.413 and -1.505 MeV (corresponding to the fragmentation schemes 2 and 3 respectively of the lj'15/2 orbital) 
since the two C(j = i 3 / 2 ) values obtained from these equations (for either C( J= i 5 / 2 ) = -1.413 or -1.505 MeV) were 
different. On the contrary, the Cy =13 / 2 ) values obtained by evaluating Eqs. (12) and (13) with Cq - =1 5 / 2 ) = 1.004 
MeV (associated to the fragmentation scheme 1 of the lj 15 / 2 orbital), are close (7.131 and 7.156 MeV). Substituting 
the average of these C(j = i 3 / 2 ) values (7.143 MeV) and C(j = i 5 / 2 ) = 1.004 MeV into Eqs. (12) and (13), we obtained 
the values C ^- 7 ) = 4.074 MeV and AC( f _ 7j = 6.139 MeV that are close to the ones specified at the beginning of this 
paragraph. These new AC”^ = 7 ) and C^~ 7 ) values are displayed as squares in Figs. 1 and 2. The discarding of negative 
C(j = i 5 / 2 ) values is consistent with the general tendency observed in the C)y values corresponding to single-particle 
neutron states of being positive. 

The GAMOW code [23] was run for the lji 3 / 2 orbital, using the GWS+SO potential with C(j = i 3 / 2 ) = 7.143 MeV . 
A positive binding energy of 4.473 MeV was obtained for this unbound orbital which represents an excitation energy 
above the 2g 9 / 2 level of 8.410 MeV in 209 Pb. In odd nuclei, such as the 209 Pb nucleus, some excited states can be 
described as admixtures of single-particle states or quasiparticle states coupled to collective excitations of the even- 
even core [26], in this case the 208 Pb nucleus. Previous works support the idea that the coupling of a single-particle 
state with surface vibrations is mainly responsible for the damping process of the single-particle mode [27-29] . Based 
on the so called quasiparticle-phonon model (QPM), Giai et al. [26] found (Fig. 4 of Ref. [26]) that the strength of 
the j'i 3 / 2 state above 7 MeV excitation energy in 209 Pb depends strongly on its coupling with the first 3“ and 5~ 
vibrational states in 208 Pb. In particular, a bump obtained at 8.8 MeV excitation energy was shown to be due to 
the coupling with the 5“ state. Similar calculations were also carried out [26] for the kn /2 and hn/ 2 states in 209 P 6 
and a QPM inclusive spectrum that included the strength functions of the k 17 / 2 , ji 3 / 2 and hn / 2 states was generated 
and compared to the experimental 208 Pb(a, 3 He) 209 Pb singles spectrum measured by Beaumel et al. [30] at 120 MeV 
incident energy and 0°. The centroids of the theoretical distribution of transfer cross sections between 6 and 12 MeV 
excitation energy in 209 Pb, were shown to be shifted to higher energies by approximately 0.5 MeV in comparison with 
the centroids observed in the experimental distribution. This indicates that the excitation energy value predicted by 
us at 8.410 MeV for the lj 13 / 2 state may be related to the theoretical 209 Pb fragment obtained by Giai et al. [26] at 
8.8 MeV or represents the average excitation energy of all the fragments calculated between 6 and 12 MeV. 

The calculation of the AC(^ =5 ) value corresponding to the open triangle displayed in Fig. 1 is based on the 207 Pb 
fragmentations (see Table III) measured by Gales et al. [21] for the \h 9 j 2 and 1 hn/ 2 orbitals. This data point is too 
far from the solid curve. On the contrary, the C^~ 5 ) value, shown also as an open triangle in Fig. 2, is consistent 
with the C( vs. £(£ + 1) value predicted by the solid curve. This agreement facilitated the calculation of new C( :j 
values for the h orbitals that are different from those displayed in Table IV. For this purpose, we substituted into 
Eqs. (12) and (13), the C(£ =5 ) value corresponding to the open triangle of Fig. 2 and the Avalue predicted 
by the solid curve displayed in Fig. 1 (open diamond). We then run the GAMOW code [23] with the new Cij values 



TABLE IV: Ct,j values of the surface potential Vsu (Eq. (11)) that make coincide the theoretical and measured binding 
energies of the shell model states analyzed in the current work. Ct and A Ct are the average and the difference respectively of 
any pair of Ctj values corresponding to a particular quantum number t. The results labeled by a, b and c are the Ctj values 
obtained for the fragmentation schemes 1, 2 and 3 respectively shown in Table II for the lj’ 15/2 orbital. 

orbital (nlj) l Ctj (MeV) Ct (MeV) A Ct (MeV) 


4Sl/2 

0 7.509 

7.509 

0 

3pi/2 

1 -3.307 

-2.842 

-0.930 

3^3/2 

-2.377 



3^3/2 

2 5.626 

5.536 

0.180 

3^5/2 

5.446 



2/5/2 

3 -4.631 

-4.522 

-0.218 

2/7/2 

-4.413 



2^7/2 

4 3.031 

2.063 

1.934 

2^9/2 

1.096 



17*9/2 

5 7.055 

1.436 

11.238 

17*11/2 

-4.183 



1 *11/2 

6 3.866 

1.562 

4.608 

1*13/2 

-0.742 



l/l5/2 

7 1.004“ 




-1.413 6 

-1.505“ 


to determine the corresponding binding and excitation energies that are shown in Table V. Although the excitation 
energy centroids E n (j displayed in this table for the I/19/2 and l/;-n/2 orbitals (3.590 and 8.543 MeV respectively) differ 
from those contained in Table III (4.014 and 8.077 MeV respectively), the new centroids lie inside the experimental 
excitation energy regions where Gales et al. [21] measured the fragments associated to the lh 9 / 2 ( E x = 3.415 to 5.620 
MeV) and l/in /2 (E x = 7.010 to 9.220 MeV) orbitals. The very strong background present in the spectrum region 
where the fragments of these two h orbitals were detected [21], suggests the convenience of carrying out new strength 
measurements to clarify the disagreements mentioned here. 

Using Eqs. (12) to (15), the strength Cpj of an orbital characterized by j and £ can be written as 


Cy — C(j = p±i/2) — + l)] 2 + /3(±)£(£ + 1) + 7(±) > (16) 

where a (±) = 0.0058 MeV, [3 (+) = -0.4422 MeV, /3 ( _) = -0.328 MeV, 7 (+) = 7.6496 MeV and 7( _) = 7.3904 MeV. 
The convenience of expression (16) is to have an easy way of determining the Cpj value required in Eq. (11) to obtain 
the correct binding energy. Eq. (16) applies to all the orbitals studied here except the p and f orbitals shown as solid 
triangles in Figs. 1 and 2. It would be advisable to revisit the investigation of these two orbitals in the future: on the 
one hand, the experimental publication [22] from which the data associated to these orbitals was extracted, seems to 
be solid and well supported [25]. On the other hand, the A Cp vs. £(£ + 1) and Cp vs. £(£ + 1) curves depicted as 
solid lines in Figs. 1 and 2, impressively reproduce the data corresponding to £ = 0, 2, 4 and 6, and make reasonable 
predictions for £ = 5 and 7. We note, without currently having a good explanation for it, that the p and f triangles 
lie close to dashed curves obtained only by shifting the solid lines shown in Figs. 1 and 2 by 1.1142 and 8.9355 MeV 
respectively. Based on this, the strength C'^ of the p and d orbitals (£ = 1 , 3) can be expressed in terms of Cpj as 

= C(j_£±i/2) = C(j—e± 1 / 2 ) — 5(±) 3 (17) 

where 5 + = 8.3784 MeV and <5_ = 9.4926 MeV. 


4. SUMMARY 


We have used a generalized Woods-Saxon + a spin-orbit potential to reproduce the experimental binding energies 
of single-particle and single-hole neutron orbitals existing in the region around the doubly-magic 208 Pb nucleus. These 
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TABLE V: Binding energies ( B n tj) of l = 5, 7 orbitals and corresponding excitation energy centroids ( E n tj ) in 207 Pb and 209 Pb 
respectively. The calculation of the Cej values is explained in the text. 

i orbital (n t j) Ctj (MeV) B n (j (MeV) E n tj (MeV) 


7 IJ 15/2 

1.004 

-2.513 

1.424 

lj’13/2 

7.143 

4.473 

8.410 

5 U 111/2 

0.046 

-15.911 

8.543 

1 ^ 9/2 

2.826 

-10.958 

3.590 


orbitals are contained in neutron shells that extend from TV = 82 to 126 (single-hole) and TV = 126 to 184 (single¬ 
particle). A remarkable agreement has been found between the binding energy values obtained by us for the 4 si/ 2, 
3d 3 / 2 , 3d 5 / 2 , 3p 1( / 2 and 3 p 3 / 2 orbitals using the WS + SO part of the GWS + SO potential only, and those determined 
in a previous investigation by Wang et al. [6] following a novel theoretical approach to solve the Schrodinger equation. 
The WS and SO parameterization and associated parameter values were taken from a recently published work [5]. 

Similarly to a common outcome reported in the literature [1-6], the WS + SO part of the GWS + SO potential did 
not produce by itself a global agreement between theoretical and experimental energy values. We have accomplished 
this task with the addition of the SU potential. For each orbital characterized by £ and j, we determined the strength 
C^j—£- ti/ 2 ) of this potential that makes the theoretical energy (calculated with the computer program GAMOW [23]) 
coincide with the experimental value (determined from neutron transfer reaction experiments carried out on 20S Pb 
[19-22]). 

A challenging result found in this investigation is the quadratic dependence of the SU potential strength x/ 2 ) 

on the quantity £{£ -(- 1). The physical origin of this dependence needs to be understood. It is puzzling that while 
the C(j = £- ti/ 2 ) values corresponding to the s, d, g, h, i and j orbitals are nicely described by the quadratic expression 
C(j=£±i/ 2 ) = a(±) [£(£ + l)] 2 + P(±)£{£ + 1) +7(±)> the strengths of the p and f orbitals (Cv^ii/ 2 )) are relatively close 
to the predictions of the shifted expression = C(j = g-ti/ 2 ) ~~ <5(±)- The different behavior exhibited by the p 

and f orbitals in comparison with the general tendency followed by most of the orbitals needs further investigation. 

It would be important to find out if the quadratic £{l+\) dependence of C(j = £.hi/ 2 ) obtained for 208 Pb, also applies 
to single-particle and single-hole neutron orbitals existing around other doubly-magic nuclei such as le O, 40 Ca, 48 Ca, 
56 Ni, 100 Sn , 132 Sn and, if this is true, determine if a global quadratic function with coefficients parameterized in 
terms of the neutron magic numbers, could be generated. It would also be worthwhile to perform similar investigations 
related to single-particle and single-hole proton orbitals. 

The explicit dependence of the SU potential strength on £(£ + 1) allowed the prediction of the binding energies of 
orbitals located inside and outside the neutron shells analyzed. The excitation energy obtained for the IU. 3/2 orbital 
may be related to the excitation energy of a lj' 13/2 fragment or the average of various fragments calculated by Giai 
et. [26] with the quasiparticle-phonon model where each excited state originates from the coupling of a single-particle 
mode with a vibrational state of the core. Based on this, we encourage the search of a formal connection, if there 
is any, between the £[£ + 1) dependence of the SU potential and the quasiparticle-phonon description or any process 
that may clarify the physics behind the SU potential included in this investigation. It is clear that further work is 
required to understand formally from a microscopic point of view the interaction mechanism responsable for the SU 
potential. 

Although the excitation energy centroids obtained for the h orbitals studied here differ from those reported previ¬ 
ously [21], their locations are inside the excitation energy region where the relevant 207 Pb fragments were detected. 
Based on the strong background present in the spectrum where the fragments were observed [21], the remeasurement 
of the fragment strengths is advisable. 

To the best of our knowledge, this may be the first time that the binding energies of single-particle and single-hole 
neutron states around 208 Pb are succesfully reproduced. This has been possible by the addition of a surface poten¬ 
tial (quadratically dependent on £(£ + 1)) to the traditional Woods-Saxon plus spin-orbit potential. In a previous 
investigation explained in the book by Nilsson et al. [31], the addition of a potential term (linearly dependent on 
£{£ + !)) to the harmonic oscillator potential make possible the reproduction of the level ordering of the finite square 
well potential. In the same work, an additional inclusion of a potential proportional to the dot product £ ■ s resulted 
in a successful reproduction of the single-neutron and single-proton level ordering measured from the spectra of nuclei 
near 208 Pb. In spite of the achievement mentioned in Ref. [31], the agreement between calculated and experimental 
energies was not satisfactory. 

It is well known that the harmonic oscillator basis has been extensively used in shell models and mean field theories 
for spherical and deformed nuclei [6]. Due to the asymptotic behavior, the use of this basis, especially in cases such 
as exotic nuclei, is questionable. The determination of the wavefunctions of a more realistic potential like it may be 
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the one used in this research would be a valuable task in future investigations. 
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